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We present a brief review of our recent efforts to develop a FDR-preserving field theory 
for the stochastic dynamic density functional model, emphasizing the essential structure of 
the theory. 

In this proceeding paper we review our recent efforts of constructing a field theory 
consistent with the fluctuation-dissipation relation (FDR) for the stochastic dynamic 
density functional model which can describe dynamics of colloidal systems (without 
hydrodynamic interaction) as well as classical atomic liquids. Such a field theory can 
provide a systematic theoretical basis for the standard MCT 1 ) which involves totally 
uncontrolled approximations such as factorization of the 4-body density correlations, 
and hence can offer a viable possibility of systematically improving the standard 
MCT. Here we hope to elucidate the essential structure of our theory. For further 
details we refer the reader to our recent papers. 2 ) 

Our starting point is the following stochastic time evolution equation for the 
density field p(r, t) 

d t p(r, t) = V ■ ( P (t, t)V-^U + r,(r, t) (0-1) 



Sp(T,t). 

where the Gaussian thermal noise rj(r, t) has zero mean and the variance 

< 7?(r, t)ri(r',t') >= -2TV • Ur, i)W(r - r')l S(t - t') (0-2) 

where T is the temperature of the system, and the Boltzmann constant ks is set to 
unity. In (|U-1|) . F[p] is the effective free energy density functional which takes the 
following form: 

F[p] = T J dr p(r) ( In ^ - l) + i J dr J dr' 8p(r) U{r - r') 5p(r') (0-3) 

where 5p(r, t) = p(r, t) — po is the density fluctuation around the equilibrium density 
Po- in (|Q~3j) the first term is the entropic contribution i^[p], and the second one the 
interaction contribution Fi nt [p\. 
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The above stochastic model (10- lj) - (jQ-3[) was proposed by one of the present au- 
thors 3 ^ as a mesoscopic kinetic equation which can describe the long-time dynamics 
of both colloidal particles and classical liquids. The equation for the coarse-grained 
density p(r,t) was obtained (with U(r) replaced by — Tc(r), c(r) being the direct 
correlation function of the system) via adiabatic elimination of the fast-decaying 
momentum field in the fluctuating hydrodynamic equations 4 ) of dense liquids. For 
this case, (j0 3f) takes the Ramakrishnan-Yussouff (RY) form. 5 ) The same form of the 
exact Langevin equation was derived 6 ^ for the microscopic density p(r, t) of Brown- 
ian particles with bare interaction potential U(r). For further discussions regarding 
the nature (and controversy) of the dynamic equation (|0-ip - (|0-3p . we refer to Ref. 7 ) 

The Fokker-Planck (FP) equation corresponding to (jO-ip and (|0 2|) can be written 

as 

where P([p],t) is the probability distribution of the density variable p(r). The equi- 
librium distribution P e [p] oc exp ( — F[p]/T) is a (possibly the only) stationary 
solution of the FP equation (I0-4D . which is rendered possible due to the prescription 
of the multiplicative noise correlation (|0 2|) . 

Note that the noninteracting part of the free energy gives rise to the linear 
diffusion for the density fluctuation 

V . f p ( r ,t)v|^4) =TV 2 p (0-5) 



Sp(r,t) 

This is an interesting fact since the two nonlinear contributions, Fid[p\ and the extra 
factor of p(r, t) in (|0-5p generate the linear diffusion which is physically expected for 
the noninteracting particles. In order to generate physically correct linear diffusion 
for the system of noninteracting particles, these two elements of nonlinearity are intri- 
cately connected. Actually, as discussed below, all the complications in the analysis 
of the original Langevin equation (IO-ip -( |(>3|) come from these two ingredients: 

• Nonpolynomial nonlinear density denpendece of F^lp] 

• Multiplicative noise structure, i.e., the extra factor of p(r, t) 

It has been explicitly shown by Miyazaki and Reichman 8 ) that due to the second 
element the response function takes an unusual form (see (|0-9fl ) . and that the FDR 
does not hold in a Martin-Siggia-Rose (MSR)-type 9 ) renormalized perturbation the- 
ory (RPT). Andreanov, Biroli, and Lefevre (ABL) 10 ) elucidated the origin of the 
incompatibility of the FDR with the usual RPT by focusing on a physical symmetry 
of the action integral, namely, the time-reversal (TR) symmetry. ABL pointed out 
that the origin of this inconsistency lies in the first ingredient, the nonpolynomial 
nonlinearity of i^[/o] entering into the TR transformation of the fields leaving the 
action invariant. ABL then proposed the introduction of the conjugate pair of aux- 
iliary fields to linearize the TR transformation, which guarantees preservation of the 
FDR order by order. The present development outlined below is a modified version 
of ABL's auxiliary field method. 

The action integral S[p,p]^ from which one can systematically obtain the cor- 
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relation and response functions of the density variable, is derived as 

S[[p,p}} ee f drj dt{ip[d t p- V- (pV^)] -Tp(Vp) 2 } (0-6) 

where the auxiliary field /5 is a real field, and the last term involving the quadratic p 
comes from the average over the multiplicative thermal noise r\. Here and elsewhere 
the functional dependence on p as a function of r is denoted as [p], but the functional 
dependence on p as a function of r and t is denoted as [[jo]]. In deriving (|0-6p . 
employing the ltd calculus makes the Jacobian of the transformation constant. (The 
dynamic action of this form with the RY free energy functional was first written 
down in. 11 ) Here we use the original form ()0-3[) for F[p]. ) 

The action (|0-6p becomes invariant when p and p transform under TR as 

p(r,-t)=p(r,t), p( r ,-t) = -p(r,t) + ^^\r (0-7) 

Tdp{r,t) 

The same action also becomes TR-invariant under another field transformation 

-*) = /?( r > *)) P( r > -*) = t) + ih(p(r, *)) 
V • (p(r, i)V%(r, t))) = ±^(r, *) (0-8) 

Note that both (|0-7p and (|0-8p are intrinsically nonlinear transformations: the first 
one is nonlinear owing to the noninteracting contribution Fid[p], whereas the second 
one due to the multiplicative nature of the original Langevin equation. As elucidated 
by ABL, 10 ) the nonlinear nature of these transformations is the underlying reason 
why the FDR, obeyed by the action, is not preserved order by order in the RPT 
developed for the action (|U-6|) . 

The FDR 12 ) is a fundamental relationship between the equilibrium fluctuations 
and the linear response to external perturbation. The response function R(r,t;r't'), 
describing how much the average density (p(r,t)) changes under application of an 
external infinitesimal field h e (r',t r ) (added to -F[/o]), is given by the following form 8 ) 

R(r, t; rV) = i(p(r, t)V ■ (p(r f , t')Vp(r', t')) ) (0-9) 



The unconventional form due to the extra factor p(r',t') in (|Q 9[) reflects the multi- 
plicative nature of the original Langevin equation (jO-ip and (jQ 2j) . The FDR then 
easily follows from the action (|0 6|) and the TR transformation (|Q7p using an identity 
(p(r,t)(5S[[p,p]}/6p(r',t'))) = 0: 

- ^d t C(r -r',t- t') = -R(r - r', *'-*) + R(r -r',t- t') (0-10) 

where C(r — r',t — t') = (5p(r,t)5p(r ! ,t')} is the density correlation function. The 
FDR (j0-10p is more directly obtained from the second transformation (|0-8p . There- 
fore the original Langevin equation is consistent with the FDR. 
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Though the original equation is compatible with the FDR, the RPT developed 
from the action (|0-6p was shown to be incompatible with the FDR. 8 ) It was subse- 
quently shown 10 ) that the intrinsic nonlinear nature of the TR transformations (|Q- 7[) 
and (|0-8p gives rise to this inconsistency. In order to resolve this inconsistency, it 
is thus crucial to develop a method which can properly linearize the nonlinear TR 
transformations. 

From now on, we focus on the first TR transformation (|0-7p . With the form of 
the free energy given in (|0-3p . the transformation (|0-7|) is explicitly written as 

p(r,-t) =p(r,t), p(r,-t) = -p(r,t)+ik*5p(r,t)+iT(5p{r,t)) (0-11) 

where K * 5p(r,t) = f dr' K(r — r')Sp(r',t) with the kernel K(r) = (5(r)/po + 
J7(r)/T). The nonlinear part T(5p{r,t)) in (|0-lip is the contribution coming from 
the higher (than quadratic) terms of F; L( i[p\: 

W ,,i)) = i(M.f|!l)..fi(. Wl , ()W . (0-12) 

where F^[p] = ^— J dr (5p(r,t)) 2 is the Gaussian part of F^p}. If one entirely 
neglects T(5p(r, t)) in (|0-lip . which is tantamount to approximating F^p] by its 
Gaussian part F$[p], then (jO-lip becomes linear and consequently the FDR would 
be preserved in the RPT order by order. However, F^[p] fails to generate solely the 
linear diffusion since 

V • (p(r, i)vf2^f) = TV 2 p(r, t) + -V ■ {SpVp) (0-13) 
V Sp(r,t)J po 

The second nonlinear term in (|U-13p gives rise to a totally spurious contribution in 
the RPT, which would yield incorrect nontrivial result even in the absence of particle 
interaction. 8 )' 10 ^' 13 ) It is therefore clear that the full logarithmic form of Fid[p] is 
required to obtain the correct behavior for the noninteracting case. 

The transformation (jO-lip can be properly linearized by introducing new vari- 
ables into the action. In particular, we introduce a new auxiliary field #(r, t) 

»(,, t )=^ W r,f))-I*^W-*M (0-14) 

T5p{r,t) po 

Note that HO- 14|) limits the new variable 9(r,t) solely to the nonlinear part of the 
transformation. This slight modification of the original ABL approach enables us to 
eliminate all the unphysical features which the ABL theory was stymied by. 

Now incorporating the new field (|0-14p . we obtain the ideal-gas contribution to 
the body force as 




+ —V • (SpVp) + poTW 2 9 + TV • (5pV9) (0-15) 
Po 
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Comparing (|0-15j) with (|0-5p . we note that the sum of the last three terms in AO- 15f) 
should vanish, and hence the new field 8(r, t) plays the role of eliminating the surpu- 
rious nonlinear term (T/po)V- (5pVp) in f)0- 13|) . This cancellation was vital to obtain 
the correct dynamic behavior for the noninteracting case nonperturbatively. 14 ) 

The resulting new action S new [[p, p, 9, 9]] can be decomposed into its Gaussian 
part Sn ew [[p, p, 9, 9]] and non-Gaussian part p, 9, 9]]: 

S new [[ P , P, 9, 0}] = Sg w [[p, p, 9, §]} + S™[\p, p, 9, §}} 
Sg w [[p,p,9,§}] = j dv j dt [ip[d t p - TV 2 p - p TV 2 9 - p Q V 2 U * bp 

-Tp (Vp) 2 + i99} 

S™[[p,p,9,9]] = J " dvj dt{ip[-V-(6 P VU*5p)--V-(6pVp) 

- TV ■ (5pV9) - T5p(Vp) 2 - i9F{8 P )\ (0-16) 

where #(r, t) variable appears from the exponentiation of the delta-functional con- 
straint 5[9(r,t) — T{p{r, t))] . The new actions Sg w [[p, p,9,9}] and <S^[[p, p, 9, 9]] 
are now separately invariant under the linear TR transformation 

p(r, -t) = p(r, t), p(r, -t) = -p(r, t) + iK * Sp(r, t) + i9(r, t) 
9(r,-t) = 9(r,t), §(r, -t) = 9(r, t) + id t p(r, t) (0-17) 

It is easy to show that the modulus of the associated transformation matrix O is 
unity (det O = — 1). Though with the constraint (|0T4p the three underlined terms 
in (|U-17p sum to vanish, each should be kept in the explicit calculation of RPT since 
their presence is crucial for separate invariance of the actions Sg w [[p, p, 9, 9]] and 
Sneu>[[Pi Pi 9, 9]] under TR, which enables one to construct the FDR-preserving RPT 
from these actions. Nontheless, as physically expected, the ultimate effect of these 
three underlined terms should be their cancellation, as explicitly demonstrated in 
the one-loop calculations in our work. Note that the linearization (|0-17p inevitably 
generates the nonpolynomial (logarithmic) nonlinearity —i9J-{5p) in the new action, 
while the original action (|0 6j) contains polynomial nonlinearities only. That is the 
price to pay for linearizing the transformation and thereby preserving the FDR in a 
RPT. 

Now the response function also takes a new form as 

R(r,t;r',t') = ^5p(r,t)9(r',t')) (0-18) 

The FDR (jOTOp then immediately follows from the linear transformation (|0T7p by 
taking correlation of the last member of (|U-17p with idp(r,t)/T. 

One is now ready to develop a FDR-preserving RPT for the new action (j0-16l) . 
A full account can be found in our recent paper. 2 ) The dynamic equations for the 
correlation and response functions are formally given by the matrix Schwinger-Dyson 
equation which involves the self-energies. These self-energies can be systematically 
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computed from the two-loop and higher-order two-particle irreducible diagrams with 
vertices mapped out of the non-Gaussian part of the action (|Q16p . Since in this work 
we are aiming to describe the dynamics of fluctuations around the equilibrium state 
(we are excluding the crystalline state from the equilbrium state), the equilibrium 
information must enter as input into the theory through the initial values of the 
dynamic correlation functions. Therefore it is crucial to provide the correct static 
input. Through this procedure, we are able to write down the dynamic equation for 
the Fourier transformed density correlation function C(k, t) alone: 



d t C(k,t) 



To(k) Jo 



S(k) 



S^(k,t-s)d s C(k,s) (0-19) 



where ro(k) = S(k) / \poTk 2 ) is the 'bare' life time, with S(k) being the static struc- 
ture factor, and E^(k,t) and E-g(\t, t) are the self-energies. Note that HQ- 19[) is 
a nonperturbative closed dynamic equation for C(k, t) only since the self-energies 
can be expressed solely in terms of C(k, t). This was demonstrated in the second 
reference of 2 ) and will be exemplified below. 

The one-loop expressions for the two self-energies in (10- 19|) are given by 



Epp(Kt) 



T 2 



^ 2 (k,q) 



k 2 
Po 



F(k,q) C(q,t)C(k-q,t), 



£ ~(M = -^ / F(k,q)C(q,i)<7(k-q,i), 
l Po J<\ 

V(k, q) = [(k • q)c(q) + (k • (k - q))c(k - q) 



(0-20) 



where c(k) = 1/po — l/SQs.) is the direct correlation function of the system. A 
salient feature of the present theory is that due to the correct static input, the bare 
interaction potential U(r) entirely cancels out, and the direct correlation function 
naturally emerges in the one- loop expression (|0-20j) . 

Another intriguing aspect of the this theory is the fact that the nonperturbative 
equation (|Q- 19|) takes a distinct form from that of the standard MCT in its structure: 
the first term in the convolution integral in (|0- 19|) does not involve the time derivative 
of C(k, t). Moreover, it is very likely due to this fact that the equation (|U-19p with 
the one-loop self-energies ()0-20p becomes unstable in the long time region. This kind 
of situation has been encountered in the mode coupling approach to the dissipative 
stochastic systems. This undesirable feature can be eliminated by noting that the 
conventional memory function can be reexpressed in terms of the so-called irreducible 
memory function. 16 )' 17 ) In particular, the projection operator approach leads to the 
exact dynamic equation for the correlation function of the chosen slow variable A(t) 



as 



dtC A (t) 



-\E A \C A (t)+ f dsM A {t-s)C A {s), 
Jo 



C^(z)=C A (0) z + \E A \-M^(z) 



(0-21) 
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where \E A \ represents a short time relaxation rate in the system, and M A (t) is the 
conventional memory function. In (|0-2ip . C A (z) etc are the Laplace transforms of 
C A (t) etc: C A (z) = J °° dt e~ zt CA(t), etc. The memory function M A (t) turns out to 
be further reducible to the irreducible memory function M A rr {t): 

M A (t) = MT(t) - \E A \~ l f dsM A {t - s)M A rr (s), 

Jo 

* rL,irr / \ 

M \{z) = A {Z . (0-22) 

l + \E A \^M L / rr {z) 

Note that M^ trr (z = 0) can grow indefinitely when the global relaxation time grows 
indefinitely in contrast to M A (z = 0). The above two eqs. lead to the dynamic eq. 
for C A (t) 

d t C A (t) = -\E A \C A (t) - \Ea\~ 1 f dsM A rr (t-s)C A (s), 

Jo 

For dissipative systems with detailed balance, while the mode coupling approxi- 
mation directly applied to the ususal memory kernel M A (t) in (|0-2ip often leads 
to unphysical results, the same approximation applied for the irreducible memory 
kernel M A rr (t) in (|0-23p always yieds physically sensible results. See. 17 ) 

Adopting the irreducible memory function approach, we rewrite (10-19j) into the 
form ^ 

d t C(k,t) = -^-C(k,t) - [ dsM(k,t-s)d s C(k,s) (0-24) 

Jo 

The irreducible memory kernel Ai(k,t) then obeys the following equation: 18 ) 

M(&t) = Eft%t)+E p6 %t)+ j dsM(k,t-s)Z pp (k,s) (0-25) 

where U pp (k, t) = U pp (k, t) / [p^Tk 2 ). Note then that the sum of the two self-energies, 
that is, the first and second terms on RHS of (I0-25P yields the standard mode coupling 
kernel: 



£ MC (Kt) = S pp (k,t)+U d§ {k,t) 



T 

Wo 



2 



(k-q) C (q) + (k-(k-q)) C (k-q) C(q, t)C(k - q, t) (0-26) 



where k = k/k. Now, when Ai(k, t) is obtained by solving (|0-25p iteratively, the 
convolution integral generates the terms ds UmcO*-^ ~ >s)^' / 3 / o(k, s) + • • • . All 
these terms belong to the higher-loop orders and can be dropped in the one loop 
order, that is, M.(k, t) — > EmcO^i t) retains the one-loop two-particle irreducible 
structure. It is therefore perfectly legitimate up to the one-loop order to retain the 
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first two terms, ignoring the terms generated by the convolution integral on RHS 
of (|0-25p . Substituting A1(k, t) = UmcQ*-, t) into (|0-24p . one recovers the standard 
MCT equation: 



We emphasize again that the bare potential completely cancel out in (|0-27p with 
(|0-26j) as well as in the previous equation (|0-19p due to the correct static input 
relation. We believe that this cancellation is no accident and can be traced back 
to our use of the Legendre-transformed vertex functional where bare interaction 
potential no longer appears explicitly. 

Here we discussed how to construct a FDR-consistent field theory for a dynamic 
density functional model applicable to colloids as well as to atomic liquids, with 
emphasis on the underlying structure of the theory. The time-reversal symmetry 
of the action integral played a vital role in developing such a field theory. The 
consistency with FDR is guaranteed when the time-reversal field transformation is 
linearized through the introduction of a new set of auxiliary field variables. The 
renormalized perturbation theory is then developed for the new action incorporating 
these new variables. Reappearance of the logarithmic nonlinearity in the new action 
is the price one has to pay for the linearization. 

For noninteracting particle systems we recover a simple diffusion law as one 
expects, which is not the case in some recent works. 8 )' 10 )' 13 ) For interacting particle 
cases, the dynamic equation for the density correlation function with one-loop self- 
energies is not only different from that of the standard MCT, but also is likely to 
be subject to the long-time instability. The irreducible memory function approach 
enables us to recover the standard MCT in the one-loop order of the renormalized 
perturbation theory. Another salient feature of our theory is the full cancellation of 
the bare interaction potential and the emergence of the direct correlation function 
which embodies the effect of the interaction potential, due to the correct static input. 

The present theory can be extended to several directions. We can first attempt 
to perform higher-order loop calculations. For this case, however, since the present 
renormalized perturbation theory is not an expansion in terms of the smallness pa- 
rameter, it is much desirable to find such a smallness parameter in order for the 
higher-loop calculations to be really meaningful. Despite some difficulties already 
encountered before, 21 ) there is a possibility one can develop such a perturbation 
theory for system with the Kac-type long-range interaction 22 - 1 employing the inverse 
force range as a smallness parameter. It would be also interesting to extend the 
theory to the inhomogeneous system 23 )' 24 ) under an external field where the spatial 
translation invariance breaks down. The calculation of the multibody (like 4-point) 
correlation functions 25 )' 26 ) can also be carried out within the present field theory 
formulation. We can also attempt to extend our approach to the genuine nonequi- 
librium problems such as ageing dynamics 27 ) and to dynamics of colloids under 
external shear flow. 28 )' 29 ) 

One profound problem related to the field theory approach to the glassy dynam- 
ics is how to incorporate the thermal activation processes into the present MSR-type 



PoTk 2 
S(k) 




(0-27) 
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field theory. Actually the dynamic density functional equation (10-lj) was proposed 3 ^ 
as a new mesoscopic kinetic equation which contains not only the nonlinear feed 
back mechanism leading to the standard MCT, but also can naturally incorporate 
the thermal activation processes. The latter statement is based on the argument 
that the only stationary solution of the FP equation (|0-4p is the equilibrium state, 
provided that the FP operator is nonsingular. 30 ) That is, the system will certainly 
remain ergodic since it will relax to the equilibrium state. Yet our one-loop analysis 
totally misses thermal activation, which indicates that thermal activation is a non- 
perturbative process. So far we have only used the uniform stationary state of the 
action in our loop expansion. We may have to consider some symmetry-breaking 
stationary configuration of the action around which a loop expansion can be made. 
It is a truly challanging problem to incorporate thermal activation process within a 
MSR-type field theory. 
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